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Recently Kleinert and Schulte-Frohlinde calculated renormalization group (RG) 
functions for the cubic model in (4 — e)-dimensions within the five-loop approximation 
of the renormalized perturbation theory |]l|. The critical (marginal) dimensionality 
Nc of the order parameter field found on the base of those expansions turned out to 
be smaller then three in three dimensions. This means that the critical behavior of 
the model should be controlled by the cubic fixed point with a specific set of critical 
exponents. Thus, calculation of the critical exponent values for the cubic model is an 
actual problem. Since the model of concern describes the critical behavior of cubic 
and tetragonal crystals undergoing magnetic and structural phase transitions, critical 
exponent estimates can be compared with experimental data. 

The field-theoretical RG series are known to be badly divergent. To extract the 
physical information relevant to predicting the critical behavior of real systems they 
should be processed by a proper resummation procedure. The series for Nc obtained in 
Ref. |jl|] proved to be alternating in signs that allowed to resumm it by the simple Pade 
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method. At the same time, the critical exponent series are irregular and cannot be 
processed by the ordinary (Pade, Pade-Borel, Pade-Borel-Leroy) techniques. A more 
sofisticated method of the Borel transformation combined with a conformal mapping, 
although regarded as the most universal procedure, is inapplicable as well, because 
it requires to know the exact values of the asymptotic parameters characterizing the 
high order behavior of the series. Nowadays those parameters have been evaluated 
for the simplest case of the 0(A^)-symmetric models only ^, and calculating them 
for anisotropic models is a most difficult problem as yet unsolved. As an exception 
one should mention the anisotropic quartic quantum oscillator which represents a one- 
dimensional (/j'^-field theory with the cubic anisotropy. For the perturbation expansion 
of the ground state energy of this system the asymptotic parameters were found in 
Ref. 1^. Within the assumption of the weak anisotropy the large-order asymptotic 
behavior of the /3-functions for the cubic model was deduced in Ref. and then used 
for determination of the stability of the cubic fixed point in three dimensions @]. 

The aims of the present work are i) to suggest a new approach to summation 
of the divergent field-theoretical series, which is based on the Borel transformation 
combined with a conformal mapping, but which does not involve the exact values of 
the asymptotic parameters. This approach, being the result of a computer study of 
the Borel summation procedure, is tested on functions expanded in their asymptotic 
power series, on caclulation of the ground state of the unharmonic oscillator, and on 
estimation of the critical exponent values for the basic models of phase transitions, ii) 
to obtain critical exponent estimates for the cubic (simplest anisotropic) model in three 
dimensions from the record five- loop e expansions using the developed technique. 

We should like to emphasize that the problem of processing divergent series arises in 
various fields of physics, where the perturbation theory is employed but the parameter 
of expansion is not small. So, developing a resummation procedure which could be 
effective where conventional methods fail is of general interest. 

A modification of the Borel procedure via a conformal mapping was introduced 
in Ref. and used for processing series J2kfkg^, whose coefficients at large order k 
behaved as k\k^°{—aQ)'' . To the series J2k fkQ^ the function 

F{g-a,h) = r e-^^(-)'d(-)B{x) (1) 

is associated. The Borel transform B{x) is the analytical continuation of its Taylor 
series X^fc a^T{b+k+i) ^^ absolutely convergent in the unity circle. Usually a = ao, whereas 
the parameter b may be related to the exact asymptotic value in a variety of ways 
[0, ||. Conformal mapping lo = transforms the cut-plain C\[— 1, —oo) onto the 

unity circle, and the semi-axis [0,cxd), the domain of integration, goes over into the 
interval [0,1). Supposing B{x) may be continued over the cut-plain, the composite 
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function B{x{lij)) is holomorphic within the unity circle and its Taylor series admits 
integration before summation. In order to eliminate possible singularity of B{x{ll!)) at 
a; = — 1 an additional parameter A was introduced in Ref. pl and used in Refs. j^, |T0|| : 
B(x{uj)) = ■ It is chosen from the condition of the most rapid convergence of 

the series 

that is from minimizing the quantity |1 — ^i{g'a\) I ' ^h^^^ L is the step of truncation 
and F^^g; a, h) is the L-partial sum of the series @). In practice, one deals with a piece 
of the series only, where the asymptotic regime might not be reached. For this reason, 
in Refs. [0, |ll|, the parameter h was varied in a neighborhood of 60 • We believe 
that in the case of the unknown exact asymptotic value ao similar manipulations may 
apply to the parameter a as well. 

Our approach to using the procedure described above consists in the following. 
While processing the asymptotic series of a priory given functions we have revealed 
that the quantity ^^((7; a, h) remains stable as the parameters a and h vary in a wide 
range. With the "number of loops" increasing, the (a, {))-dependence becomes weaker 
and weaker, and the difference between Fi^g^a^h) and the exact value F{g) can be 
made whatever small. This observation enables us to employ the Borel transformation 
with a conformal mapping for processing series whose exact asymptotic behavior is 
unknown. Thus, we postulate stability of the result of the processing with respect to 
variation of a and 6 as a basic principle underlying our approach. 

Let us demonstrate how this principle works on some simple examples. At first, 
consider a model function 



00 

,.4 



k=o 



whose coefficients fk behave as ^ ^ at large k. For each a and b in (0) we find the 



k 

set of values {A^j„(a,6)}, at which the quantity |1 — -^^"^f^f^l^ as a function of A, 
reaches its local minimums. To each element of that set corresponds the particular 
value J-'l{g; a,b). It can be shown that taking into account relative contributions of 
lower order terms of the series {k < L) when optimizing the A-parameter just weakly 
affects the final result, therefore it is sufficient to minimize the relative contribution of 
the last accounted term only. At Fig.l (a) a few curves J-'l^{0.5; ao,b) are presented, 
which, as calculations show, deviate from the exact number S{g) = 1.48631082... by 
no more than 5 ■ 10~^ within the range < 6 < 18. 
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For the function J-'{g) the exact asymptotic value of a is ao = 4. As the parameter 
a shifts from the picture does not change quahtatively, ahhough at large numbers 
one can see a reduction of the 6-stability interval, as displayed at Fig.l (a). For small 
a the set {A^j„(a, fe)} is just empty within the interval [0,6]. In all cases the lines in 
their stability domains are at the same level with the mentioned above accuracy. So, 
on the graph of J-'l^{g; a,b) depending on two parameters a and b there would be a 
horizontal part - plateau rapidly destroyed at its boundaries. 

The size of the plateau depends essentially on g. When g gets smaller the stability 
domain is enlarged. That is in accord with intuitive expectations that the decrease of 
the expansion variable should result in better convergence. So, the calculations show 

does 



.FJ^(0.1;ao,fe)-J^(0.1) 



that for the family of five curves ^{3(0.!; ao, b) the relative error , ^-^^ -^-^ 
not exceed 2 ■ 10"^ within the interval < 6 < 25 (see the table in Appendix). On the 
contrary, the growth of g leads to reduction of the stability domain. Let us follow how 
its size changes from loop to loop, assuming g = 1. Presented at Fig.l (b) are the most 
stable curves J-'l{1.0; Oq, b) for L = 6, 9, 12. As seen from the graphs, with accounting 
more terms of the expansion the behavior of the curves improves, and the dispersion 
of the values goes down. The graph illustrating convergence of numerical estimates of 
J-'{g) for = 1 depending on the truncation number L is depicted at Fig. 2 

When processing functions {xdx)^° j^dx ~ Yj'kLo{~^Yk\k^° g^ , where bo > 

0, we observed the similar stability with respect to a and b. We also processed the 
six-loop pieces of RG series for the critical exponents of the 0(A^)-symmetric model in 
three dimensions and found very weak dependence of the output on the transformation 
parameters. The numerical values of the critical exponents computed proved to be in 
a good agreement with those of Ref. [0. 

Using the proposed summation method for processing the series of the ground state 



energy S{g) of the unharmonic oscillator [jT3[ with the Hamiltonian H = + gx'^ we 



observed the same behavior of the corresponding curves as for the model functions 
considered above. In Table 1 we present the estimates of S{g) at g = 1 depending on 
the length of the series. 

TABLE I. Numerical estimates for the unharmonic oscillator ground state energy 
at g = 1. 



L 


8 


9 


10 


11 


12 


Exact value 


^(1) 


1.392376 


1.392357 


1.392344 


1.392349 


1.392351 


1.392352 



Note that for L = 8 our estimate is by one order closer to the exact value than the 
number 1.391655 ± 0.004562 found in Ref. on the base of Wynn's e-algorithm. 

The numerical analysis fulfilled allows to apply the summation method introduced 
to finding numerical estimates of the critical exponents for the (4— e)-dimensional cubic 
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model [|1| in three dimensions (e = 1). The 6-dependence of the results of processing 
the exponent 7"^ for = 2 at various fixed a is presented at Fig. 3. The parameter a 
ranges from 0.2 to 1.5 with step 0.1. Distinct oscillations correspond to the small values 
of a and b. As a grows, the behavior of the curves becomes smoother and the extended 
horizontal interval appears reaching its maximum length at about a = 0.5. Further 
increase of a makes the stability interval getting shorter and causes the rapid growth 
of the curves at its boundaries. It is essential that for all values of a the horizontal 
parts of the curves are at the same level. Averaging the results of the processing over 
the stability domain gives the number which we adopt as 7"^ sought for. The accuracy 
for this approximate value is determined through the dispersion due to the variation 
of a and b. Similar behavior of the curves was observed for other critical exponents. 
Applying to them the developed algorithm we obtain the critical exponent values for 
the cubic model for various listed in Table 2. 

TABLE 11. Critical exponents for the cubic model in three dimensions from the 
five-loop approximation in e for some A^. 



N 


V 


V 


7 


7sc 


7s c 


2 


0.0350 ±0.0003 


0.6277 ±0.0010 


1.2358 ±0.0040 


1.2334 


0.19 % 


3 


0.0375 ±0.0005 


0.6997 ±0.0024 


1.3746 ±0.0020 


1.3732 


0.10 % 


4 


0.0365 ±0.0005 


0.7225 ±0.0022 


1.4208 ±0.0030 


1.4186 


0.15 % 


5 


0.0358 ±0.0004 


0.7290 ±0.0016 


1.4305 ±0.0040 


1.4319 


0.10 % 


6 


0.0354 ±0.0003 


0.7301 ±0.0016 


1.4322 ±0.0040 


1.4344 


0.15 % 


00 


0.0350 ±0.0003 


0.7108 ±0.0010 


1.3993 ±0.0020 


1.3967 


0.19 % 



The case N = 00 corresponds to the Ising model with equilibrium magnetic impurities. 
In this limit the Ising critical exponents a and i/ are renormalized according to Fisher 
[|T^. The values of r], u, and 7 were computed by processing original series for r], z^~^, 
and 7~^, while 7sc was found by the scaling: 'jsc = i^(2 — 1]). It is seen from the table 
that 7sc differs from 7 by no more than 0.2 %, and this may testify the good accuracy 
of the estimates obtained. 

Let us compare our results with earlier estimates of the critical exponents for the 
pure Ising model. It is known that due to the symmetry of the initial Hamiltonian 
of the cubic model at A^ = 2 the critical exponents for the cubic and the Ising fixed 



points coincide. There were found the following estimates for the Ising model 
r] = 0.035 ± 0.002, 1/ = 0.628 ± 0.001. These numbers are in excellent agreement with 
the data of Table 2 at A^ = 2. Our estimates are also in a good accordance with the 
results of Ref . |]12| , where a substantially different method was employed. 

In conclusion let us formulate the results of the present paper. A new approach 
to summation of divergent series has been suggested. The method employes the Borel 
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transformation combined with a conformal mapping. It relies upon the stabihty of the 
result of processing on the transformation parameters and therefore does not require 
to know the exact asymptotic behavior of the series. The method has been tested on 
the functions expanded in their asymptotic series and applied to estimating critical 
exponent values for the cubic model. The principal observation is that within the 
new approach to summation of the perturbative series of both simple and complex 
(anisotropic) models exhibit the same behavior. This allows one to apply the devel- 
oped technique to process divergent series arising in a number of anisotropic models 



describing phase transitions in real substances |]T6[. It can be expected that the new 
summation method may be useful in other fields of physics (e.g. QCD and QED) 
where one deals with divergent series and conventional resummation techniques are 
inapplicable. 
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Figure 1. a) Graphs of dependence of ^2(0.5; a, 6), L — 13, on parameter h for some 
values of parameter a. 

b) Graphs of dependence of T1{g]l,h) on parameter h in sixth, ninth, and twelfth 
orders in g for g — 1. 

Figure 2. Convergence of the numerical estimates of the function 

^— oo 

depending on the approximation order L for g = 1. The dashed line corresponds to the 
exact value 1.368427... The upper and lower broken hues yield, respectively, the upper 
and lower boundaries for the estimate. The best estimate is given by the middle line. 

Figure 3. Curves demonstrating dependence of the result of processing the critical 
exponent 7"^ on parameter b for various values of a for the cubic model at A/^ = 2 in 
five-loop approximation. 



A Appendix 

Table of the results of processing the asymptotic series for the function 

/+00 2 . 
-00 

with the following values of the parameters: g = 0.1, a = = 4, L = 13, X E [—2,4]. 
Numbers J^l^ig) correspond to different X''{g; a, b) realizing the most rapid convergence 
of the processed series. The exact value is JF(O.I) — 1.674085856.. . 
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